PERRON THEOREM IN THE MONOTONE ITERATION METHOD FOR 
TRAVELING WAVES IN DELAYED REACTION-DIFFUSION EQUATIONS 



AMIN BOUMENIR AND NGUYEN VAN MINH 

Abstract. In this paper we revisit the existence of traveling waves for delayed reaction diffusion 
equations by the monotone iteration method. We show that Perron Theorem on existence of 
bounded solution provides a rigorous and constructive framework to find traveling wave solutions 
of reaction diffusion systems with time delay. The method is tried out on two classical examples 
with delay: the predator-prey and Belousov-Zhabotinskii models. 



1. Introduction 

We shall be concerned with the existence of traveling waves solutions for the delayed reaction 
diffusion system 

(1) *iM = jD fHM + /(ut ). 

at ax z 

Due to their important applications in population dynamics and biological models, see e.g. [271 [5J 
ll[^ll[3[ll[M[Sl[li[Ml[31[2ni[Iil[Il[3Il[37], equations such as © have evolved from the 
simple one dimensional scalar reaction diffusion equation to systems that include a delay in time 
for more realistic modeling. Traveling waves, although are a classical topic in applied mathematics, 
remain a driving force in the study of (TTJ) . We refer the reader to [101 130j which contain many 
surveys on methods used to study traveling waves in parabolic differential equations. Note that 
when delay is introduced, most methods would fail if they are not modified appropriately. 

Recently, Wu and Zou [32] have extended the method of monotone iterations to deal with 
the existence of traveling wave fronts for delayed equations ([lj (see also [20l [21] [14l [15]). This 
produces a monotone sequence of positive functions that converges to a traveling wave front, which 
is an increasing positive solution of an equation of the form 

(2) D0"(O-c0'(C)+/c(^)=O, eei. 

This method is contingent on the construction of a pair of upper and lower solutions, which satisfy 
the inequality versions of ([2]). It turns out that this issue is closely related to the existence and 
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uniqueness of smooth and bounded solutions of nonhomogeneous equations of the form 

(3) Dx"(t) - cx'(t) -0x(t) +g(t) = 0, for all t G R, 

where the constant c, /3 > 0. In differential equations results of this type are known as Perron 
Theorem, see [TJ [T31 [531 123 and for interesting applications see [5J [TJ and the references therein. 

We would like to point out that the smoothness required for upper and lower solutions of 
([2]) obstruct the search for bounded solutions of ([3]). One faces the following dilemma, see e.g. 
36, 32, 20J . An excessive relaxation of the smoothness of upper and lower solutions of ([2|) simplifies 
their finding, but would not generate the sought monotone iteration scheme. The main reason is 
due to the failure of Perron Theorem for this class of weak solutions of ^ . 

In the light of the above remark, the question of existence of traveling waves front in predator- 
prey or Belousov- Zhabotinskii models with delay, as addressed in [20], remains open. 

The purpose of this paper is two folds: First to set up a rigorous framework for the Monotone 
Iteration Method and then apply it to the predator-prey and Belousov- Zhabotinskii models with 
delay. 

We now briefly outline the plan of this paper. In the next section, to recall the main concepts 
and tools, we discuss a modified version of Perron Theorem for C 1 -solutions of Next we show 
how crucial it is to the monotone iteration method. Remarks and counter-examples are used to 
explain the pitfalls of non-smooth upper solutions. In Theorem [TT] one finds a rigorous framework 
to construct fail-safe upper and lower solutions which are then used for the delayed predator-prey 
and Belousov-Zhabotinskii equations, see Theorems fT5| [T71 

To conclude we would like to emphasise that Perron Thereom dictates C 1 -smoothness which 
makes the search for upper and lower solutions much harder than C°-smooth solutions as in |20j , 
and of course easier than C 2 solutions, which seems to be impossible for the above mentioned 
equations. 

2. The Monotone Iteration Scheme 

In this section we introduce a modified version of Perron Theorem which is central to a rigorous 
framework for the monotone iteration scheme. We show that this framework provides a clean 
procedure for constructing upper and lower solutions of reaction diffusion equation with delay. 

2.1. Bounded solutions of nonhomogeneous equations. The theory of bounded solutions 
of nonhomogeneous equations is a classical topic of the theory of ordinary differential equations, 
which can be found in [TJ [13l [T9l [25] and the references therein. 
Let us first consider the equation 

(4) u"{t) + au'{t) + (3u(t) + f(t) = 0, t e M, u(t) € K, 

where / is a function that is continuous and bounded on R\{0} and has the right and left limits 
at x = 0, / (0 + ) and /(0~); We always assume that a and /3 are real numbers with f3 < so that 
the characteristic equation 

A 2 + aX + 13 = 
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has two distinct real roots of opposite signs Ai < and A2 > 0. To our knowledge, all available 
results on classical solutions deal with bounded and continuous forcing terms / on the entire real 
line (see e.g. [TJ |T3l [19l [25] ). As we allow jump discontinuities in / in {U, we need to modify the 
concept of solutions as well as the conditions for their existence and boundedness. We first need a 
definition. 

Definition 1. Suppose that f is a bounded and continuous function on R\{0} and both /(0 + ) and 
/(CP) exist. Then, a function u defined on R is said to be a generalized solution of if 

(1) u and u' are bounded and continuous on R, 

(2) u" exists and is continuous on R\{0}, and both u"(0~) and it"(0 + ) exist. 
Below is a version of Perron Theorem for generalized solutions with discontinuous /. 

Lemma 2. Consider equation with [3 < 0, and assume that 

(1) f is a bounded and continuous function on R\{0} and both f(0 + ) and f(0~) exist, 

(2) Equation holds in the classical sense for all t except possibly at t = 0. 
Then, Eq. (^J) has a unique generalized solution u given by 

(5) u{t) = Gf{t) := ^-L_ e x ^f(s)ds + e x ^f(s)ds^ . 

Proof. Observe that function Gf, which is defined by ([5]), exists and is bounded. A simple compu- 
tation shows that Gf is in fact continuously differentiable on R, and (Gf)" exists and is continuous 
on the whole R with a possible exception at t = at which both (Gf)"(0 + ) and (Gf)"(0~) exist. 
It is important to observe that while Gf and (Gf)' are bounded on R, Gf is a particular solution 
of (H]) on each of the two disjoint intervals (—00, 0) and (0, +00). Obviously, v — u — Gf is then a 
classical solution of the homogeneous equation associated with (|4]) on (0, 00), and hence 
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(t) - Gf(t) = ae Xlt + be X2t for t> 



where a, b are constants. From the boundedness of u and Gf on (0, 00) it follows that b = i.e. 

u(t) - Gf(t) = ae Xlt for t > 0. 
Similarly the classical solution u — Gf on (—00, 0) is of the form 

u(t) - Gf(t) = de X2t for t < 0. 
Combining both behaviors we have 

ae Xlt , t > 0, 
de X2t , t < 0. 



(6) u(t) - Gf(t) 



Use the fact that although u(t) — Gf has no second derivative at t — 0, it is still continuously 
differentiable on the real line, and in particular at t = 0, to deduce the interface conditions 

{a = d 
Al<2 = A2C?. 

Since 7^ Ai 7^ A2 7^ 0, we must have a = d = 0, that is u — Gf = 0. This completes the proof of 
the lemma. □ 
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Remark 3. Lemma\^ shows that the continuity of u' plays a crucial role in the uniqueness of the 
solution. The failure of the continuity even at a single point, would lead to a non-uniqueness of 
bounded solutions. We now illustrate this important fact by a simple and yet explicit example. 

Example 4. Consider bounded solutions for the equation 

(7) y"(t) - y(t) = where t € E. 

Since f = 0, by Lemma\Mthe only bounded continuously differentiable solution (that is, generalized 
solution) would be 0. However, if we allow for a non- differentiable function to be a solution, then 
the continuous function 



(8) y{t) 



t < 0, 



satisfies ^ for all t £ ]R\{0}. Obviously, it is a nonzero bounded and continuous function on K, 
but not differentiable at t = 0. 

Remark 5. The function y in |3J) can also serve as a counter-example to the identity (see the 
proof of [20l Lemma 3.5],) 

t />oo 



(9) ___ ( / e ^ t - s >cp(s)ds+ I e^ {t - s > V {s) 



1 rn 

V(t) + 1 — T E e^- T >\y'{T ]+ )-y'{T-)) 



A2 — Ai 

^2 e Mt-T s ) iyl{T . +) _ yf{T _ ))t 



j=k+i 

k 



3=1 

where y is a piece-wise C 2 solution of 

dy"{t) - cy\t) - (3y{t) = <p(t) on R\{T , ...,T m }. 

In the particular case when d = 1, Ai = — 1, A2 = 1, c = 0, (3 = 1, m = 0, Tq = 0, the expression 
(Q) defines a solution to y"(t) — y(t) — for all t 7^ 0. Since (p — 0, the left hand side of (0) is zero 
while the right hand side is not 

_j m k 

= 2/(«)+2 E e- t (y'(0 + )-2/ , (0-))--E et (y'(0 + )-y'(0-)) 

3=fe+l 3=1 
k 



1 I " 

= y(*)+2 E e- t (-2)--Ee'(-2) 
j=fc+i ]=i 

= y(t)- e - t +e l ^0 



We hope that that the argument for a continuous y' is by now clear and self-evident. 



PERRON THEOREM AND MONOTONE ITERATION METHOD 



5 



2.2. Delayed reaction diffusion equations. Consider the following system of reaction-diffusion 
equations with time delay 

an\ du{x,t) d 2 u{x,t) 

(10) __=.D__ +/(„,), 

where t > 0, x £ R, u(x,t) £ R", D = diag(d%, ■■■,d n ) with dj > for all i — 1, „., n, / : 
C([— t, 0], R") — > R" is a continuous, and Ut(a;) is an element of C([—r, 0],R"), defined as 

u t (x) = u(t + 9, x), 9 £ [-T, 0], t > 0, x £ R. 

Throughout this paper we will assume that 

(If) /(6) = f(K) = and f(u) ^ 0, for u € R" with < it < K, 

where 6 denotes the function ip £ C[-t, 0], R") such that ip(6) = := {0, 0, 0} T € R™ for all 6 e 
[-t,0] and denotes the function ^ e C[-r,0],R n ) such that if>(9) = K := (K U K 2 , -,K n ) T E 
R n , with given positive Ki, i = 1, 2, n.. In order to use comparison arguments, we use the natural 
partial order in R" to compare two vectors x — (xi, ...,x n ) T , y — (yi, ...,y n ) T £ R", that is x > y 
if and only if Xi > yi, for all i = 1, n, and if there is at least an i in {1, n} such that Xi < yi, 
we write x < y. The "interval" [0, K] consists of all vectors v £ R" such that < Vi < Ki for all 
i = 1, 2, n. 

As usual we look for traveling wave solutions of (|f 0|) in the form u{x,t) — <j>(x + ct), where 
G C 2 (R, R"), and c > is a constant. Substituting u(x,t) — cj>{x + ct) into (fT0|l leads to the 
following wave equation 

(12) £>0"(e)-^'(O+/c(^) = o 1 eeR, 

where f c : X c := C([-cr, 0], R") R™, defined as 

/ c W=/(n il> B (9):=^(cB), 9£[~r,0}. 

Throughout this section we assume that there exists a matrix (3 — diag{(3\, ...,/?„) with (3 > such 
that 

(13) /cW)-/cW+^(0)-#)]>0, for all Wel c , < ip(8) < <f>(6) < K, 

where := (0, 0) T £ R" and K := {K, K) T £ R". 

The main purpose of this paper is to look for solutions 4> of (|f 2p in the following subset of 
C(R,R") 

r := {ip £ C(R,R") : (p is nondecreasing, and lim (p(£) = 0, lim </?(£) = K} 

^ — > — oo £ — ^+oo 

The solution (j> is called a monotone wave front of (fTT))) . Next we define an operator H : BC(R, R") — -> 
BC(R, R") by 

(14) H(<p)(t) = f c ((j> t ) +(3</>(t), 0eC(R,R n ). 

The following monotonicity lemma was proved in 32J. 



Lemma 6. Assume that \13\) and \1J$ hold. Then, for any <f> £ T, we have that 

(1) H(<f>)(t) > 0, t £ R, 

(2) H(cj))(t) is nondecreasing in t £ R, 
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(3) H(<ip)(t) < H(4>)(t) for all t G R, iftp G C(R,R") is given so that < ip(t) < <j){t) < K for 
all t G R. 

Our definition of upper solutions, which is found below, requires more smoothness and bound- 
edness than those in [32] and [20] . 

Definition 7. A function p G C 2 (R, R") with p, p' , p" being bounded on R is said to be an upper 
solution of if it satisfies the following 

(15) Dp"(t) - cp'(t) + f c (p t ) < 0, for all t G R. 

Definition 8. A function p G C (R, R") is said to be a quasi-upper solution of U2\) if 
(1) 

sup||p(t)|| < oo, sup ||p' (f) || < oo, 

teR teR 

(2) p"(t) exists and is continuous on R\{0}, and 

sup \\p"(t)\\ < oo, 



tei 

(3) lim t ^ - anc ^ lim ( _ +0 + /°"(^) ea; * s ^ 

(4) p{i) satisfies 

(16) D/(t) - cp'(i) + / c ( ft ) < 0, /or a// t G R\{0}. 

Similarly, we define the concept of lower and quasi-lower solutions to (|12p by switching the 
inequalities in (JTSJ) and (16]). Let G C(R,R n ) such that < (f>(t) < K. By Lemma H we have 
< H(cj))(t) < K, so H((j)) is a bounded function. Note that 



2d~ <°' A * : = 2d~ >0 

are the real roots of the equation 

di\ 2 - cX - Pi = 0, i = l, 2, n. 

The existence of a unique generalized bounded solution of 

(17) Dx"(t)-cx'(t)- f3x{t)+H(cf)(t) =0, for all teR 
follows by Perron Theorem, see Lemma [2] is 

G(H(<P)) = {G 1 {H 1 {<j>), G 2 (ff 2 (0)),..., G n (H n (J>))) T 

where 

(18) Gi{H(4>)) := - 1 . ( f e^-^Hii^ds + / + °° e x ^-^ H^^ds) . 

Ah — A 2 i \J-oo / 

The Monotone Iteration Scheme is constructed as follows: We start out with a quasi-upper 
solution </>o, and then use the recurrence formula 

(19) <f> n := G(H (</> n _i)), n=l,2,.... 

Lemma 9. Assume that (f>,<po £ T are, respectively, a quasi-lower and quasi-upper solutions of 
{H), smc/i that 4>{t) < 4> (t) for all t G R. Tften 
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(i) 0i G r for all t G M, 

(ii) 0i is a upper solution of flip , and 

<f>(t) < <t>i{t) < <fo{t), for all tern.. 
Proof. For the proof of (i) we refer to J32[ Lemma 3.3]. Next, we prove (ii). Clearly <f>\ = 

G(H((j> )) G C 2 (K,R n ) and by Lemma[H <j>\ is also the unique bounded (classical) solution of the 
equation 

Dx"(t) ~ cx'(t) - fjx(t) + H(<p Q )(t) = 0, for all t G M. 

By Lemma [51 H{4>o) > 0, and so 

D0'/(i) - c0i(t) - < for all t G K, 

which means that <f>\ is a upper solution of (|12p and from its definition <f>i > 0. We now show that 
the sequence {4> n } is decreasing. To this end, set 

w{t) :=0i(t)-£)(i), 

(20) r(t) := -Dw"{t) - cw'(t) + 0w(t), t G M. 

By (|19[) and the assumption that 0o is a quasi-upper solution of (|12p we see that r,(t) is non- 
positive and bounded for every i = 1, ...,n and £ G R\{0}. Moreover, u)j(-) and the i-component 
of (|20| satisfies all conditions of Lemma O Therefore, it follows that as a bounded generalized 
solution of ([20|) . Wi satisfies 

Wi(t) = Gr t {t) := (J e^-^r^ds + e x ^-^n(s)ds^ < 

which yields w(t) — (f>\ — 4>q < 0, i.e. <pi < </>q- Next we show that <fi < 4>\. To this end, we again 
set 

v(t) 

(21) s(t) := -Dv"(t) - cu'(i) + /3w(t) > 0, (el. 
By the definition of quasi- lower solution and <pi and by Lemma [2] we have 

(22) = G (*) (t) := (|* e^Ct-O^^dg + e^-^R) > 

which implies <pi(t) > <f>(i). □ □ 

Remark 10. As shown in Example^ without assumption on the continuity of the derivative of 
</>0; formula \22i) may not be true because we know only that w is bounded and it may not be of C 2 . 



Below we will use the notation Bq 0jK ](l,R") := {g G 5C(K,R n )| < g(t) < K, Vt G R} as a 
closed convex subspace of the Banach space BC(R, M") equipped with the sup-norm. We are now 
ready to state the Monotone Iteration Method. 

Theorem 11. Assume that 

(1) there exist a quasi-upper solution of ilty) 4> Q G T and a non-zero non- decreasing quasi-lower 
solution 0o such that < <fio(t) < <f>o(t) for all t G K, 



8 



AMIN BOUMENIR. AND NGUYEN VAN MINH 



(2) H is continuous on BC[ .k] (R, K"), and 
(23) sup ||i%)|| < oo. 

Then, the following assertions hold: 

(i) The sequence {4>n}^Li, defined as above, is a decreasing sequence in T, 

(ii) linin^oo 4> n {t) = 4>(t) is a monotone wave front of Moreover, this limit is uniform on 
each compact interval of the real line. 

Proof. By Lemma [9] Claim (i) is clear. To prove (ii) we use Arzela-Ascoli Theorem. Note that by 
Lemma El all {(f> n }%Li are C 2 functions. From 

-( f e x ^ t ^HM(s)ds+ f +0 ° e^-^HM^ds) , 

2i \J-oo Jt J 



we deduce 



(24) \{4>n + i)'S)\ < 



2sup v3Gr || J ffM|| 



|Aij — A2i| 

Therefore, on each interval [-N, N], the set {4> n }%Li is equicontinuous and by Arzela-Ascoli 
Theorem, we have a uniformly convergent subsequence {<^n fc }fe on [-N, N]. Since the sequence is 
monotone we can conclude that the sequence {4> n }^Li converges uniformly to a continuous and 
non-decreasing function <f> such that < <p(t) < K on every compact interval of the real line. Since 
(f> € 5C[o,K](K,t") and H is continuous on BC[ 0i k](R,R"), for every fixed (t,s), we have 

lim e^-^HMis) = e^-^H^is). 

Lebesgue's Dominated Convergence Theorem, then yields 

= lim (<t> n+ i(t))i = - 1 - ( f Hi(<l>)(8)da + [ °° e^-^Hii^ds) 

= (G(H(cj>))(t))i- 

This shows in particular that tj> is a classical solution of the equation 

D<t>"{t) - c<t>'{t) - (3c/)(t) + H(<P)(t) = 0. 
It remains to see that <f> £ T. First since 4>(t) < <fi n (t) for all n and t E M. 

< limsup0(i) < lira sup <f) n (t) = 0. 

t — > — oo t — > — oo 

Next to show that lim^oo <fi(t) = K, we use ^ 4>o(t) < 4>(t) for all t, and the fact that 4> is 
non-decreasing means that limt^oo <f)(t) exists as a vector L € [0,K]. Obviously, L ^ 0, and so 
< L. We now use the delay and as in [32l Proposition 2.1], we must have / C (L) = 0. But the 
assumption on the function / makes it impossible if L ^ K. So, lim^oo (f>(t) — K, and 4> € T. □ 
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3. Traveling Waves for the Predator-Prey Model 
Consider the predator-prey model with diffusion and delay in time 



(25) 



= d^^A +ru(x,t) [(1 - ±u{x,t)) - av(x,t)] 
^at'^ = d 2 d a^'^ + v(x, t) \-v + bu(x, t-r)), 



where u{x, t), v(x, t) are scalar functions, [x, t) G K x [r, oo), and di,d 2 , r, P, a, b, v, r are all positive 
constants. Define R- valued functions fi(4>, ip) and f2{4>, tp) as follows: For each (j>, ip G C([— r, 0],R), 
the functionals 



fi(<f>,i>) = r(j){0) 



1 



0(0) 
P 



oV>(0) 



/ 2 (0,V)=^(Q) (-i/ + 60(-t)). 



If we assume that 



1 > 



<^ P > 



P6 b : 
then the model has a non trivial positive steady state 

T 



v 1 



(1 



Pb 



The above model with r = has been treated in [27, 4, 5 (for related results see [22, 23, 24 .) 
The case with delay r > was first considered in [20] . As noted in Remark [5] the smoothness 
was overlooked in the proof of [20, Lemma 3.5]. This gap is due to the failure of Perron Theorem 
for this type of "suppersolutions" . Hence, the question on the existence of traveling waves for this 
model with delay remains open. A quick remedy is provided by the results of the previous section, 
which require us to construct quasi-upper and quasi-lower solutions which have more smoothness 
than "suppersolutions" in [20] . 



3.1. Wave equations and upper, lower solutions. Let us agree to denote by (fx(x + ct) 
u(x,t), ip 2 {x + ct) :— v(x, t), and £ = x + ct . We recast ([23]) into 



C ^(0 = dM'(0 + rMO [(i - - <m(Q 

cy 2 (0 = da^'fO + pa(0 [-z/ + b Vl (Z - ct)] = 0, 







p 



a<^ 2 (£) 



which can be re- written as 

(26) di^(0-<vi(0 + m(0 

d2^ / (O-c^a(O + Va(O[-" + ^i«-c7-)]=0- 
We will find monotone solutions of the above wave equation such that 

(27) 

(28) lim ip 2 (t) = 



lim ifi(t) = lim ipi(t) =p 

t — > — oo t — >+oo 



lim tpn(t) = q 

t — >+oo 



V 

b 
1 



v \ 

„ \ ~ ~Pb) 
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We now check the quasi-monotonicity of (/i,/2) T , that is, there are positive numbers 0i,02 
such that 

/i(^i> 2 ) - fiWi, *h) + Mm - ^(o)) > o 

/ 2 (^i> 2 ) - Mipi,^h) + P(<h(o) - Mo)) > o 

for all 

0<^)<^)<K:= (^(l-^y 

Using the definition of the functions (fi,f 2 ) T we can easily come up with a pair of positive numbers 
/3i , [3 2 that satisfies 

/ /3i>max(-r + ^ + f,0) 
02 > f. 

3.2. Upper solutions. Let us define ip(t) :— ((pi (t), ip 2 (t)) T , where 



(20) Mt) = < ^ eAlt ' w*"° and ^ 2 (t) = J 2 eA ; 



where the parameters are chosen so that 

c + \/c 2 — 4dir c v 1 / v \ 

Al = 2^ ' A2= 2^' P =b> q= a\ l -pb)- 

Therefore, Ai and A2 are simply the roots of the quadratics 

(30) diA 2 -cA + r = 0, d 2 X 2 - cA = 0. 

Lemma 12. Assume that v < b, and d 2 > 2d\. Then, there exists a constant c* — c*(a, b, r, v, P) > 
such that if c > c* , (f>o(t) := (ipi(t), f 2 (t)) T is a monotone quasi-upper solution of the wave 
equation. 

Proof. It is easily seen that both <pi and ip 2 are continuously differentiable, monotone and obviously 
lim <fi(t) = lim ifi(t) = p := — 

t^ — 00 t^+oc 

lim <p 2 (t)=0 lim <p 2 (t) = q := - (l - -^r) . 

t— >— 00 a \ Pbs 

Next, we have 

^-e A i*, t < „, , f ^e Al *, t < 



2 



... f ^e A2 *, t<0 .... [ ^e A2t , t<0 



S^-A* i>() r 2 W-| ^| e -A 2tj t>Q • 

A crucial property of the functions ipi, ip 2 is their smoothness. We now look for conditions on the 
parameters so that they form upper solutions. In doing so, we examine the cases when t < and 
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t > separately. Substituting the above expressions into the first equation of (|2l)|) and using ([50)1 
yields for all t < 

pA? _ pAx p Alt _ rp _ arp? Alt A2t 

2 2 2 2P 2 

_ rP Ait o n>g Ai« Aat 

2P 2 

< 0. 

Similarly for £ > 0, we have 



"M/' -' "J ('/ - |e _A2(t_CT) 



(,-Ie-)G 



(rp - ^ - ar-pg) + g + ^ - C pAi) e~ Alt + 



arpg e _A 2t arpg e _ ( A 1+ A 2 )t 



2 4 

To check the sign of the above expression, we first factor out rp, which is positive, 
(3!) (l - I - „) + (± + f - 4) + f - 3 e -0. + M, 

A sufficient condition on P, can be obtained from (|3T|) and use if A2 > Ai, which holds if 0I2 > 2d\ 

P \b 

to be non-positive for all t > 0, is 

(32) V - < 1 

(33) \- + 2aq< 2c— 

P r 



(i 


+ a<? 


2Ar 

r 




1 

p + 


2ag - 




e~ Alt < 



For (|33|) to hold it is sufficient to take c large enough. More precisely, if c > \fAd\r, then Ai > 2c 



which, in turn, leads to c > \\ {j> + 2aq) r. Thus the value of c* can be estimated by 



1 / 1 



c* = max [ -y + 2a<jJ r, y / 4dr> 

Corollary 13. In case f < P < \, then there c* = 0. 

Proof. It is enough to see that (|31[) can also be written as 

1\ (\ Ai\ \ + /e- A ^ + e- Alt e -(A 1 +A 3 )t 
1 - t; + 777^ - c — e + W 



P ) \2P r ) * V 2 4 

1 Ai Af /e- Alt -l\ / e -A2t + e -A lt (A 1+ A 2 )^ 



2P r \ 2P 

It turns out that we have an upper solution, provided c > and 

(35) 1 -^<°> ^h P< \ 

no restrictions on di, d.2- □ 
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Now we check the second equation of the wave equation. Starting again with t < and substi- 
tuting the expressions for (pi, (p 2 into the second equation, together with (|30|) yields 



' 2 

for all i < since t — ct < 0. 



r 



V 2b 



(t-cr) 



-^e A2t (2-e Al ^ CT )) <0, 



Next, for t > 0, substituting into the second equation of the wave equation yields 
- d 2 §e~^ - c^e-^ - „ ( ? - f ) + b (g - fe"^) (p - | e -M«-cr)) 
= -(„-,*)§ + („- pb) I (e- A2t - 1) - 2cA 2(?e - A2t - P^ e -Mt-cr) + P^ e -x l{ t-c T ) e -x,t^ 
Since 

pb = v 

the above expression reduces to 

- 2cX 2 qe- x ^ - E^ e -^-cr) + P^_ e -x l{ t-cr) e -X 2 t 
< -2c\ 2 qe-^ t - ^ e -^(t-cr) + qV e - M t-cr) e -X 2 t 



-X 2 t _ ^L p -Xl{t-CT) 

2 



1 -r 



-x 2 t 



= —2cq\ 2 e 
< 0. 

and therefore the lemma is proved. 

3.3. Lower solutions. We will construct a quasi-lower solutions to (f2"o]) as follows: Set 



□ 



which is a root of 



v\ := — . 



d\v x — cv\ = 0. 



Let us define functions ip 2 (t) = for all i 6 R, and 
(36) Vi(*) = 



ake Vlt : t < 
k-ake~ Ult , t> 0, 



where the constants fc, and a are positive, sufficiently small and will be determined at the end. A 
simple computation yields 



V4W 



v 1 ake' Jl \ t < 
v 1 ake-' Jl \ t > 0. 

and the first equation of (|26[) with t < gives 

(diakvi — cakv\)e Vlt + rake Vlt 
1 



vfake Ult , t < 
~vlake~ Vlt , t > 0. 



1 - -ake Ult 



= rah 



> 0, 



1 - -^ake Vlt 



provided < ak < P. 



PERRON THEOREM AND MONOTONE ITERATION METHOD 



13 



For t > 0, substituting the derivatives of ipi into the first equation of (f2T))) leads to 
- divfake~ Ult - cu%ake~ Ult + r(k - ake~ Vlt ) 



1- -^(k-ake- Ult ) 





+ ak I 


M) 







2rk 
~P 



> 0. 



provided that c and k are fixed, < k < P, and a is sufficiently small. 
Next, we need to choose a, k so that 

^l(t) < M*)- 

In fact, a simple computation shows that we may choose k < v/{4b) and a < (4Pb)/v. 

Obviously, this function (V'i,V'2) T is a quasi-lower solution to (|26[) . Using the same procedure 
as for upper solutions we will obtain a smooth lower solution as desired. 

Therefore, we have proved 

Lemma 14. Let all assumptions of Lemma \12\ hold. Then, for every fixed c > c* , where c* is 
as in Lemma \T2\ there are (sufficiently small) constants a and k such that the function 4>o(t) := 
(i/ji(t),0) T , where ipi(t) * s defined by S36\) is a quasi-lower solution of H26\) that satisfies: 

(37) o < Mt) < Mt), t e r, 

where (j)Q is the quasi-upper solution that is defined by i29\) and mentioned in Lemma [7^1 

Finally, by Theorem [TT] we have 

Theorem 15. Assume that v < b, and di > 1d\ . Then, there exists a constant c* — c*(a, b, r, i>, P) > 
such that if c> c* , Eq. i25\) has a wave front solution u{x,f) = tpi(x + ct), v(x, t) = if2{x + ct). 
Moreover, if ¥ < P < i, c* may be chosen to be 0. 

4. Belousov- Zhabotinskii Equations 

In this section we will apply the results in Section 2 to prove the existence of traveling waves to 
Belousov- Zhabotinskii Equations 



(38) 



jtu(x,t) = -§^u{x,t) + u(x,t)[l - u(x,t) - rv(x,t - r)]; 
j%v{x, t) = -^2v(x, t) - bu(x, t)v(x, t), 

where r > 0, b > are constants, u and v are scalar functions. We refer the reader to |27j for 
more details on the history as well as applications of this kind of equations. Traveling waves in 
these models without delay were considered in various papers (see e.g. [TTJ [TBI ISi]). In the 
recent papers [5U1 [35] , the models with delay were first studied as an application of the Monotone 
Iteration Method, which was extended by Wu and Zou. As we have noted earlier, in the previous 
model, the Monotone Iteration Method reduces to finding a pair of upper and lower solutions to 
the corresponding wave equation. Again due to the failure of Perron Theorem for the concepts of 
" upper solutions" and " supersolutions" in [351 [5D] , the question of existence of traveling waves in 
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these models remained open. Below, Perron Theorem and Theorem 1 1 1 1 will be used to answer this 
question by constructing suitable quasi-upper and quasi-lower solutions to the wave equation. 

The wave equation associated with Belousov- Zhabotinskii Equations can be modified (see also 
[321 [20]) so that it takes the form 



(39) 



<p'{(t) - ap[(t) + <pi(t) ((1 - r) - ip^t) + rcp 2 (t - or)) = 
ip'Ht) - ap' 2 {t) + tyi(t) (1 - <p 2 (t)) = 0. 



We seek monotone solutions (ipi,ip2) T such that 

lim tpi(t) = 0, lim ifi(t) = 1 

t — > — oo t — > + oo 

lim tp2{t) = 0, lim ip 2 (t) = 1- 

t — > — oo t — > + oo 

To this end we recast the wave equations in the following form 

(40) yi{t) - ap'^t) + Vl (t)] - r<pi(t)(l - Mt - cr)) - p?(t) = 

y 2 '{t) - ap' 2 {t)] + b Vl (t) ~ bMVMt) = 0. 

Define the numbers Ai and /ii as 



c + Vc 2 - 4 c + Vc 2 - 46 
Ai- 2 , Ml- ^ , 

which are the roots of the characteristic equations 

A 2 - cA + 1 = 0, [L 2 - c/i + b = 

respectively. Observe that 

(41) b < 1 => Ai < /i! 

Let us define functions y>i and y>2 as follows: 



m , |e Alt , t<0, J ie^*, t<0, 



observe that Q < ifi(t) < 1 and similarly for < (fi2(t) < 1. 

Lemma 16. There exists a positive number c* = c*(b, r) such that if c > c* , then (ipi, f 2 ) T defined 
as above is a smooth monotone upper solution of the wave equation. 

Proof. First, it is easily seen that 



{^J-p Alt t < f ^ip A i* / < 

2 e ! <- _s u > ii ( ,\ I 2 ' — > 



V±e^\ t<0, .... ^-e^, t<0, 
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The second derivative exists almost everywhere, substituting the above expressions into the first 
equation of PO]) we have for £ < 

&{{t) - c^i(t) + V i(t)\ - np x {t){l -<(>2(t- or)) - <o?(t) 
= -re Alt (1 - 99 2 (* - cr)) - ie 2Alt < 

since (1 — ip 2 {t — cr)) > 0. 

For the second component we have for £ < 

^'(£) - c<4(£) + kpi® (1 - - ~ c^e^ + (\ - 



= Alt _ „\L\t\ 



a (Ai+/ii)« 



Thus by 



Ai < mi => ^'W - c</ 2 (£) + (£) (1 - <fii(t)) < -- < 



On the other hand if £ > 0, then we have 

[*>?(«) - c^(£) + ^(£)] - r^(£)(l - ^ 2 (£- cr)) - <^ 2 (£) 



e Alt - c^ie- Alt + 1 - V Alt 



- /■ ( 1 - -< A: ' 
2 



(l-^ 2 (£-cr))-^(£) 



(l - ±e" Alt ) (1 - Mt - cr)) - Xfe-^ - \e-^ 



< 0. 



Similarly the second equation yields 

^{t)-c^ 2 {t)+b Vl {t){l-ip 2 {t)) 

= zA e -^* _ c^e-^ + b{l- -e~ xj 
2 2 



bipi(t)ip 2 (t) 



-mi* + 1 (e-^t - e - A i') + -e^^ + -e~ Alt - - e ~ (Al+fll)t 



(6-c) 



-Ait 



-fJ-it 



so, when £ > 0, 

6 - c < =*> </j 2 '(£) - c^ 2 (£) + 6^1 (*) (1 - <P2{t)) < 
To summarize the situation we have 
b < c ) 



Ai < /xi J 
and by (|4Tj) the condition reduces to 



which is satisfied if 



<4(£) - c^ 2 (£) + 6v?i(£) (1 - ^ 2 (£)) < V£ G 



b < c 
b<l 



b<c 
Ai < ^1 



b < min {1, c} . 
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□ 

Finally, since the procedure of constructing quasi-lower solutions is similar for the predator-prey 
models we leave it to the reader. To conclude this section we have 

Theorem 17. There exists a positive number c* = c*{b,r) such that if c > c* , then Eq. i38\) has 
a wave front solution u(x, t) = ip\{x + ct), v(x,i) = ip2(x + ct). 
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